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The spin-resolved electronic energy band spectra, as well as partial and total density of electronic states of 
the crystal AIN, doped with Cr, Mn and Fe, have been evaluated within the projector augmented waves (PAW) 
approach by means of the ABINIT code. The Hartree-Fock exchange for correlated electrons is used to describe 
the correlated orbitals in the PAW framework. The calculated one-electron energies for electrons of spin up and 
down are very different. We have found that all the considered crystals are ferromagnetic. 
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1. Introduction 

Semiconductors doped with d-transition elements draw considerable attention of researchers due to 
their possible applications in spintronics U|. The discovery of ferromagnetism therein led to an inten- 
sification of experimental and theoretical studies. The impurities of transition elements lead to a fun- 
damental change in the electronic structure of semiconductors y]]. Crystal AIN with an admixture of 
3d-transition elements was investigated using the ultrasoft (USPP) pseudopotential approach within the 
LSDA theory 1 3] and in the LSDA+ U approximation [4]. The partial and total densities of electronic states 
of the Mn-doped crystal AIN and co-doped with B, Si, and P, with spin up and down |5], have been re- 
cently evaluated by means of the LSDA KKR program 1 6] . The parameter of the Coulomb energy U makes 
it possible to take into account strong correlations between the d-electrons of the transition element, and 
thus to improve the values of the electronic energy band spectrum of the crystal with impurity. 

However, it is known that the Coulomb energy U depends on the system. This means that it cannot 
be used in calculations in the same crystal, for example, with vacancies \7]. As it is inherent to LSDA+ U 
approach, in 01 there was also found the existence of several minima of the total energy functional 
depending on the matrix of the initial occupation of the correlated orbitals. 

The spin-resolved electronic energy band spectra, as well as partial and total density of electronic 
states of the crystal AIN, doped with Cr, Mn and Fe, have been evaluated within the projector augmented 
wave (PAW) approach |8, 9], implemented in the ABINIT 1 10] code. The exact exchange, for 3d-electrons 
only, is used to improve the on-site correlations II ill . 



2. Calculation 

The electronic structure calculations have been carried out using the PAW 1 9] method, in which the 
electron wave function with its full nodal structure \if/ n {r)) is represented in terms of the smooth nodeless 
function \xj/ n (r)) as follows: 

\VnM) = T\Vn(X)>, (2.D 
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where t is the transformation operator, 

T = i+IO02>-i#M<02i. ( 2 - 2 ) 

a a 

and \(p%{t)) is the all-electron basis function, |0£(r)) is the pseudopotential basis function, and lp„(r)) 
is projector function, which we calculated using the atompaw code 1 12]. The subscript a denotes the 
quantum numbers, and the superscript a represents the atomic augmentation sphere index. 

We have generated the PAW functions for the following valence basis states: 2s 2 2p 6 3s 2 3p 1 for Al, 
2s 2 2p 3 for N, 3s 2 3p 6 3d 5 4s 1 4p° for Cr, 3s 2 3p 6 3d 5 4s 2 4p° for Mn, and 3s 2 3p 6 3d 6 4s 2 4p° for Fe. The inclu- 
sion in the basis of 3s and 3p states of the Cr, Mn and Fe atoms slows down the calculation, but increases 
the precision and eliminates the problem of ghost states which can appear among the solutions of a secu- 
lar equation constructed on a minimal set of functions 11311 . The radii of the augmentation spheres are 1.6, 
1.3, 1.9, 1.9, 1.9 a.u. for Al, N, Cr, Mn, Fe, respectively. Substituting (IH . (2~D to the Kohn-Sham equation 



H\V>ak) = Wak)£ak (2.3) 

we obtain the system of linear equations 

T + Hr\y ak ) = T + x|i/> ak )c ak , (2.4) 

where en is a band number, k denotes the vector in the first Brillouin zone, and t + defines Hermitian 
conjugate of t. The electron density in the PAW method is determined by the sum of three terms, 

p(r) = p(r) + £[p fl (r)-p fl (r)]. (2.5) 

a 

The first term is a smooth pseudo-density, which is the Fourier series, 

P(r) = Zfnk\Vak(r)\ 2 = -^Xp(G)e iGr , (2.6) 

n,k 12 G 

where f n k is occupancy, weighted by the fractional Brillouin zone sampling volume. The terms of the den- 
sity of electrons in the atomic spheres are determined from the coefficients of the projected population 
of the states 

W ap =Lfnk<Vak\Pa)<Pp\V'ak>. (2.7) 
n,k 

Hence, the contributions of the density in the atomic sphere can be written as 

p fl (r) = X KpVi OrtyjsM. P" W = E <6<Pa W^W. (2.8) 
ap a.B 

Exchange-correlation potential was calculated in the form of PBEO [14] according to which the exchange- 
correlation energy 

£ x p c BE0 [p] = £™ E [p] + ~ (EfVseil ~E™ E [p sel ]) , (2.9) 

corresponds to the PBE exchange-correlation functional [IS], and y/ se \, p se \ represent the wave function 
and electron density of the selected electrons, respectively (lg]. The latter are the 3d-electrons of Cr, Mn 
and Fe. 

The electronic energy bands and DOS have been evaluated by means of the ABINIT code 1 10]. Integra- 
tion over the Brillouin zone was performed on the Monkhorst-Pack 1 17] spatial grid of 4 x 6 x 6. The lattice 
constant of the crystal A1N a = 4.36 A, and the supercell parameters containing 16 atoms are 2x1x1. The 
iterations were performed to ensure the calculation of the total energy of the crystal with an accuracy 
of 10~ 8 Ha. Relativistic effects are taken into account within the scalar relativistic approximation. The 
density of electronic states was evaluated from the equation 

^(£) = £<5(£-£ ak )|P;^(k)| 2 , (2.10) 
a,k 
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KpW = (Pp\Vak) = f drpf (r-R a )^ ak (r), (2.11) 
where ^^(k) are the expansion coefficients for both 

l^ak>=LOk)l<^>, (2.12) 

P 

P 

smooth and true functions, respectively, within the augmentation sphere. As can be seen from equa- 
tion <2.9t . the exchange-correlation energy functional depends on the density of electrons and on the se- 
lected wave function of strongly correlated electrons. First, the pseudo-wave function is calculated from 
equations <2.4) . Then, from equation H2.1> , the electron wave function is calculated and the density of 
electrons (2.5) is obtained. Based on the latter we find the exchange-correlation energy functional (2.9) . 
The symmetry of the crystal is described by space group P-42m (number 111) and the Bravais lattice is tP 
(primitive tetragonal). 



3. Electronic properties 

Figure [l] shows the partial and total density of electronic states of the crystal Al 7 CriN8 for a spin 
down. We see that the top of the valence band is formed predominantly with hybridized s-states of Cr 
and p-states of N. The dispersion curves shown in figure [3] indicate that for spin down, our crystal is a 
semiconductor with a direct gap. It is substantially less than the gap T-X in the usual AIN crystal. Alu- 
minum nitride in the zinc-blende structure is an insulator with a wide indirect energy bandgap of 5.3 eV. 
Such a narrowing is explained by the fact that the crystal CrN is a Mott-Hubbard-type insulator with a 
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Figure 1. The spin down partial, in l/(eV-atom-spin), 
and total densities of states for Al 7 CriN8. 



Figure 2. The spin up partial, in l/(eV-atom-spin), and 
total densities of states for Al 7 CriN8. 
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small to negligible indirect band gap [18]. The bottom of the valence band is now at the point R, while in 
the the usual crystal A1N, it is at the point F. The Fermi level is located in the band gap, quite close to the 
bottom of the conduction band. 

For spin up, we have a completely different picture, represented in figures [2] and [4] The crystal 
AlyCriNs is a semiconductor with indirect narrow gaps T-R and T-A. The Fermi level is now even closer 
to the bottom of the conduction band. The top of the valence band is formed by hybridized d-states of Cr, 
p-states of N and Al. 

We turn to the analysis of the results of calculation for the Al 7 MniNe crystal, represented in fig- 
ures[sHHl 

Figures[5]and[7]describe the electronic density of states and energy bands for spin down. From figure[5] 
we see that the top of the valence band forms the s-states of Mn and the p-states of N. Figure[7]shows that 
the crystal is a direct band gap semiconductor, with the bottom of the valence band at the point R. 

Figures[6]and[8]describe the electronic density of states and energy bands for spin up. Figure[6]shows 
that the top of the valence band is formed by the p- and d-states of Mn and the p-states of N, but the 
bottom of the conduction band forms the s-states of N. Figure[8]shows that the crystal is a semiconductor 
with indirect gap T-ZR. The Fermi level is at the top of the valence band. The presence of the spin up 
3d states of Mn at Fermi level (figure [8) and of the spin down ones in the conduction band (figured) 
qualitatively agrees with the results obtained within the DFT + U approach for Mn doped A1N fl]. 

Let us consider the results of the calculation for the AlyFeiNs crystal, presented in figures I91IT21 The 
calculation results for the spin down are presented in figures[9land[TT1 Figure [9] shows that the top of the 
valence band is formed by p-states of N, and the bottom of the conduction band consists predominantly 
of rf-states of Fe. Figure [TT1 shows the direct band gap at the T point. The Fermi level is located near the 
bottom of the conduction band. The bottom of the valence band is at the point R. 

The calculation results for the spin up are presented in figures [TOl and [121 In figure [101 it is seen that 
the top of the valence band is formed by d-states of Fe and p-states of N. The bottom of the conduction 
band forms the p-states of Al and s-states of N. From figure Il2l we see that the crystal is a semiconductor 
with indirect gaps F-R and F-A. The bottom of the valence band is at the point R. 
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Figure 5. The spin down partial, in l/(eV-atom-spin), 
and total densities of states for AlyMniNs. 



Figure 6. The spin up partial, in l/(eV-atom-spin), and 
total densities of states for AlyMniNs. 





Figure 7. The spin down band structure of Figure 8. The spin up band structure of Al 7 MniN8. 
AlyMniNg. 
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Figure 9. The spin down partial, in l/(eV-atom-spin), 
and total densities of states for Al 7 FeiN8. 



Figure 10. The spin up partial, in l/(eV-atom-spin), 
and total densities of states for Al 7 FeiNg. 





Figure 11. The spin down band structure of Figure 12. The spin up band structure of AlyFeiNs. 
Al 7 FeiN 8 . 
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4. Conclusions 

For the spin down, all the crystals are direct band gap semiconductors, with direct gap at the point 
T. And for the spin up all the crystals are characterized by indirect band gap. All the crystals show the 
ferromagnetic ordering. For AlCrN, our result is confirmed by experimental data on ferromagnetism 
detected above the room temperature Hal . For AlMnN we observe a half-metallic behavior in the sense 
that the Fermi level state density is finite for the majority spin, and zero for the minority spin, as it is 
seen in figures [7] and [8] Figure [6] shows that at the Fermi level, there are the spin-polarized carriers. 
This is well matched with the results obtained in the work 1 20] within the LSDA approach. However, our 
calculation shows that the AlCrN crystal does not possess the half-metallicity predicted in the work 1 20] . 
We have calculated the electron energy spectrum and the density of electronic states in the PAW GGA 
formalism without accounting the strong correlations of 3d-electrons. It was found that there is no state 
of half-metallicity in the AlCrN crystal which does exist in the crystals of AlMnN and AlFeN. Meanwhile, 
for the latter crystal, figures ITT1 and Il2l obtained within the PAW PBEO approach, the half-metallic state 
is not observed. As we see in table [l] the values of both direct and indirect gaps of the crystals are quite 
different. The largest contributions to the density of states at the Fermi level is provided by the states of 
d- and p-symmetry of manganese and p-states of nitrogen. The difference between the results obtained 
herein with the account of strong electron correlations, in the PBEO approach, and in LSDA theory, is due 
to inadequacy of the latter in the description of the wave functions and energy levels of 3rf-electrons. 

Table 1. The calculated electronic band gaps (for spin down, dn, for spin up, up, in eV) and spin magnetic 
moments (m in ^g) of crystals. 
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As a result, the calculation from first principles of electronic and magnetic properties of the ATTN 
systems, where T is a transition element has been made taking into account strong correlations of 3d- 
electrons. A significant narrowing of the interband gap was found for all systems. The calculated electron 
energy spectrum and the density of electronic states, as well as the value of the Fermi energy are impor- 
tant parameters for the systems considered to be essential candidates for high-temperature electronics 
and optoelectronic applications 12 ill . 

Another important factor is the issue of maintaining the ferromagnetic state at higher temperatures 
as much as possible. However, this is the subject of a separate study. 
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EiieKTpoHHi B/iacTMBOCTi Kpi/irrasia AIN, 
sieroBaHoro Cr, Mn Ta Fe 

C.B. Ci/ipoTK)K, B.M. LLlBefl 

HaLjiOHa/ibHwfi yHiBepcmeT "ZlbBiBCbKa no/iiTexHiKa", By/i. C. EaHflepn 12, 79013 m. /lbBiB, YKpaiHa 

3a/ie>KHi Bifl cniHa e/ieKTpoHHi eHepre™4H0T cneKTpn, a TaKOx napuia/ibHi vt noBHi Lui/ibHOCTi e/ieKTpoHHux 
CTaHiB Kpucra/ia AIN, /leroBaHoro Cr, Mn Ta Fe, 6y/in OTpuMaHi 3a MeTOflOM npoeKUMHMX napLjia/ibHwx XBH/ib 
(PAW) 3a flonoMOTOK) nporpawiw ABINIT. 06MiHHwR noTei-mia/i XapTpi-<t>OKa pj\n CKOpe/ibOBai-inx e/ieKTpoi-tiB 
BUKopucTOByeTbcn pj\ft onucy CKOpe/ibOBai-inx ap&nanevi b paMKax PAW. Po3paxoBaHi OflHoe/ieKTpom-ii eHeprii 
Pflft e/ieKTpoHiB 3i cniHOM Bropy i bhw3 e Ay>K e pi3Hi. Mn buhbw/iw, mo Bci po3r/iflHyTi KpwcTa/iw e c|>epoMarHe- 

THK3MH. 

lOnoMOBi cnoBa: po3paxyHOK eneKTpoHHO'i crpyKTypn, cmnbHO CKOpenbOBaHi e/ieKrpoHn, tohmavi oSmihajiR 
CKOpe/ibOBaHnx e/ieKTpoHiB, MarHimi HaniBnpoBiflHMKH, Meiop, npoeKU,ii/iHi/ix napujanbHwx XBwnb 
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